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Abstract 

We study the single-particle spectral properties of a model for coexisting 
AFM and ICDW critical fluctuations coupled to electrons, which naturally 
arises in the context of the stripe-quantum-critical-point scenario for high-T c 
superconducting materials. Within a perturbative approach, we show that 
the on-shell inverse scattering time deviates from the normal Fermi-liquid 
behavior near the points of the Fermi surface connected by the characteristic 
wave- vectors of the critical fluctuations (hot spots). The anomalous behavior 
is stronger when the hot spots are located near singular points of the electronic 
spectrum. 

The violations to the normal Fermi-liquid behavior are associated with the 
transfer of spectral weight from the quasi-particle peak to incoherent shadow 
peaks, which produces an enhancement of incoherent spectral weight near the 
Fermi level. 

We use our results to discuss recent ARPES experiments on Bi2212 near 
optimal doping. 

PACS: 71.10.Hf, 71.45.Lr, 74.25.Jb 
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I. INTRODUCTION 



The nature of the single-particle excitations in a metal close to a Quantum Critical Point 
(QCP) has become a subject of interest after the proposal that the anomalous properties of 
the metallic phase of high-T c superconductors may be related to the presence of a singular 
effective scattering amplitude, such as the one arising near a QCP of some kind of instability, 
in the Random Phase Approximation (RPA). Two realization of this scenario have been 
proposed, the first associated with an antiferromagnetic (AFM) instability ||T]-0], the second 
with an incommensurate charge- density- wave (ICDW) instability M. 

The presence of an AFM phase at zero and very low doping and of strong AFM fluctu- 
ations, revealed in neutron-scattering experiments at higher doping ||, has suggested that 
the relevant physics of high-T c superconducting materials is dominated by the AFM QCP. 
However in such scenario the QCP itself is quite far from the point of optimal doping where 
both the highest T c and the maximum violations of the Fermi-liquid (FL) behavior in the 
metallic phase are observed. To explain this evident inconsistency strong vertex corrections 
have to be advocated in order to suppress the effect of AFM fluctuations |J, which would 
otherwise be strongest at lower doping, close to the QCP. Moreover the peculiar role of the 
optimal-doping point in all classes of high-T c materials is left unexplained. 

To avoid this problem, and considering also that previous theoretical results indicate the 
presence of a charge instability, the existence of a new QCP was proposed, which controls 
the physics of the cuprates near optimal doping. Indeed, models for strongly correlated 
electrons with short-range interactions are commonly characterized by an instability with 
respect to phase separation J7|, which is turned into an instability with respect to ICDW 
when long-range Coulomb forces are taken into account @J|. 

The most direct evidence for a new QCP at (or near to) optimal doping is provided by 
resistivity measurements |J, which reveal an insulator-to-metal transition when the super- 
conducting state is suppressed by a strong pulsed magnetic field. Neutron scattering exper- 
iments evidenced charge-driven order in related compounds [frO| . Coexistence of ICDW and 
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electron gas has been detected in joint EXAFS fTT| , |T2|j and X-ray diffraction [|1^] experiments 
(see also Ref. |b§). 

The new QCP is not incompatible with an AFM QCP. The two coexist and determine 
the physics of high-T c superconducting materials in different regions of the phase diagram. 
Near optimal doping it is very likely that charge degrees of freedom play a major role. 
Moreover the dynamical charge segregation into hole-rich and hole-poor regions, associated 
with critical charge fluctuations, would enhance AFM fluctuations far away from the AFM 
QCP, due to the natural tendency of hole-poor regions towards antiferromagnetism. In this 
way a consistent framework is achieved to understand the role of spin degrees of freedom 
in optimally doped and overdoped materials. AFM fluctuations enslaved by the onset of 
charge dynamical fluctuations were evidenced in related compounds |L0 . 



Different approaches have been used so far to capture the relevant features of the single- 
particle excitation spectra in a metal close to a QCP. In the mean-field theory of spin- 
density-wave (SDW) antiferromagnetism the quasi-particle spectrum may be obtained by 
coupling the electrons to spin fluctuations characterized by a factorized susceptibility x — 
i7iS(S + l)5(uj)5(q — Q), where Q = (n/a, ir/a) is the wave-vector of the AFM structure 
and 5 is the spin quantum number of the fluctuating spin. Within this approach the gap in 
the quasi-particle spectrum is proportional to the effective coupling g. Kampf and Schrieffer 
proposed that the resulting self-energy T,(k, e) = g 2 [{s — ^-q) 1 — iiisgn(e)6(e — ^-q)], could 
be slightly modified, near the AFM QCP, due to the interaction between free electrons and 
AFM fluctuations |T5|JT6[| , which leads to a broadening |ImS| ~ g 2 T /[(s — ik-o) 2 + T 2 ] so 
that ReS ~ g 2 (e — ^-q) / [(^ — (,k-Q) 2 + ^ 2 } has a quasi-polar structure [|17]]. Thus the inverse 
scattering time at the Fermi energy is finite 1/r ~ 2r/[£|_g + T 2 } and new (quasi-)poles 
appear in the electron Green function at energies e ~ £&_q, which are usually called shadow 
bands. 

Within the framework of the nearly- AFM FL (NAFL) [[J, Chubukov et al. ||, inves- 
tigated the single-particle properties of a model for electrons coupled to spin fluctuations 
characterized by a phenomenological susceptibility of the form proposed by Millis, Monien 



and Pines |18| . An undamped susceptibility x{Qi u ) was use d to describe mainly the under- 
doped region. Within perturbation theory a topological transition of the Fermi surface (FS) 
was found at a finite value of the electron-spin fluctuations coupling constant |J, leading 
to the formation of hole pockets located around the points (±7r/2a, ±7r/2a) of the Brillouin 
zone. Such a transition is associated with the appearance of new poles in the electron Green 
function which are the result of the evolution of the electron self-energy towards a mean-field- 



like behavior similar to that suggested in Ref. flgjlql . Thus, in this scenario, the precursors 



of the AFM phase show up in the anomalous metallic phase, at higher doping, with a change 
in the nature of the quasi-particles, signalling the incipient antiferromagnetism. 

However the experimental results provide an increasing evidence for the existence of 



stripe phases with modulations incommensurate with the lattice | 19| . This makes at least 
questionable a description in terms of AFM fluctuations only. Moreover, the topology of the 
FS of optimally doped Bi2212 was recently determined with great accuracy pPf , revealing 
new relevant features, and namely an asymmetric suppression of the spectral weight near 
the Fermi level around the M points of the Brillouin zone, at points of the FS connected 
by an incommensurate wave- vector Q c = (0.47r/a, — OAir/a), in addition to a symmetric 
suppression at the same points, associated with the wavevector Q s = (it /a, ir/a). 

This may be interpreted within the ICDW QCP scenario §, as a signature of quasi- 
critical charge fluctuations close to a stripe phase. In the following we want to check the 
consistency of this interpretation, and suggest that ARPES experiments can be used to 
extract the properties of the charge-fluctuation spectrum, from the related features appearing 
in the single-particle excitation spectra. 

Since we expect that the onset of a stripe phase reintroduces AFM fluctuations, we 
consider as a starting point an effective model for coexisting charge and spin fluctuations 
coupled to conduction electrons with bare band dispersion The properties of such charge 
and spin fluctuations are generically described within a RPA around the QCP at optimal 
doping, which can be considered as the point governing the onset of the stripe phase. The 
explicit form of the parameters appearing in the resulting charge and spin susceptibilities, 



which depend on the underlying microscopic model, will be deduced below by means of a 
phenomenological analysis, mainly based on ARPES experiments. In particular, the modu- 
lation of the density waves associated with charge and spin fluctuations is characterized by 
the wave-vectors Q c and Q s respectively 

Within a perturbative approach we show that such a model exhibits a violation of the 
normal FL behavior associated with the transfer of spectral weight from the quasi-particle 
peak at an energy e ~ to dispersing incoherent bands located at e ~ 6e-Q cs , without the 
appearance of new poles in the electron Green function. The resulting bands have thus the 
aspect, but not the coherent nature of the shadow bands found within the approaches of 
Refs. P,|T5ll. Such bands could be related to corresponding broad features seen in ARPES 
experiments, such as the peak at an energy e ~ —200 meV, observed near the M points of 
the Brillouin zone in Bi2212 |21[ . 



We point out that the shadow features of incoherent nature were also found by Benne- 



mann et al. within the FLEX approximation for the Hubbard model ||22j| . The interpretation 
of the shadow bands as spectral- weight anomalies 2^] was also suggested in ref. More 
recently similar results were found within the NAFL scenario [2S|, by performing a full 
re-summation of the entire perturbation series in the high-temperature (static) limit for 
the spin-fluctuation spectrum. We point out that these results agree with our simpler per- 
turbative results, in the case of electrons coupled to spin fluctuations only. Moreover, by 
considering also the effect of charge fluctuations, we reproduce new features, which have a 
counterpart in recent ARPES experiments, and call for further investigation. 

The scheme of the paper is the following. In Sec. II we introduce the model for tight- 
binding electrons coupled to charge and spin fluctuations. In Sec. Ill we discuss the violation 
of the FL behavior for the on-shell inverse scattering time. In Sec. IV we analyze the single- 
particle spectral properties of our model. Sec. V is devoted to the comparison between our 
results and ARPES experiments on Bi2212 near optimal doping. Concluding remarks are 
found in Sec. VI. 
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II. THE MODEL 



We consider a model for tight-binding electrons coupled to charge and spin fluctuations, 
described by the Hamiltonian 

fe,cr i kq,a/3 

where = —2t[cos(k x a) + cos(fc y a)] — At' cos(k x a) cos(k y a) — //is the dispersion law for free 
electrons on a square bi-dimensional lattice for the CuC>2 planes with hopping terms up to 
next-to-nearest neighbors included, a is the lattice spacing, and /i is the bare chemical po- 
tential. In the following, for the sake of definiteness, and aiming in the end for a comparison 
with ARPES experiments, we assume that —1 < It' jt < 0, to fix the main properties of the 
bare band structure and the shape of the FS, as appropriate to the Bi2212 samples. The 
correction 5 fx ~ 0(g 2 ) to the chemical potential is treated perturbatively and is fixed to 
ensure that the number of particles is the same as that in the free-electron system with bare 
chemical potential /i. 

The second term in (P describes the interaction of electrons with charge (i = 0) and spin 
(i = 1,2,3) fluctuating fields S l , which is characterized by the coupling constants gi. The 
spin structure of the generalized electron density coupled to the fluctuating field S l _ q is given 
by the corresponding Pauli matrix r\ We assume that the phenomenological properties of 
the fluctuating fields S l _ are completely described by correlation functions of the form 

x«(g,fc0= o ,y*. , (2) 

where, for each j, Aj is some constant, flj(q) = Mj + afT + Vjjq-Qj with j q = 2 — cos(q x a) — 
cos(q y a), Qj is the characteristic wave- vector of the critical fluctuations, Uj is the inverse 
of the characteristic time scale of the fluctuating mode and the distance from criticality 
is measured by the "mass term" Mj + acjT. In particular Mj is the distance from the 
QCP at zero temperature, which we assume to depend on the hole doping per unit cell x 
and to vanish at a critical value x c , and T is the temperature in energetic units, so that 
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a, is a dimensionless constant. We point out that the cos-like form of the dispersion in 
7(g) is adopted to reproduce a q 2 behavior at small momenta, while preserving the lattice 
periodicity [Q. 

We further assume complete spin isotropy, so that all the above parameters are the same 
for j = 1, 2, 3. Thus in the following we turn to the notation j = c,s for all the constants 
related to charge and spin fluctuating fields respectively. We also introduce the dimensionless 
coupling constants A c = g 2 c A c jtv c and A s = g 2 A s /tu s . 

In the case of AFM fluctuations the susceptibility (||) corresponds to the phenomenolog- 
ical expression obtained by Millis, Monien and Pines |18| in the limit of strong damping; in 



the case of ICDW a similar expression was found close to the instability in the Hubbard- 
Holstein model with long-range Coulomb forces, within a slave-boson approach @||. In the 
present phenomenological approach the fact that charge instability leads to an enslaved spin 
modulation at the stripe QCP, is modelled by the requirement that the vanishing of the mass 
M s is guided by the vanishing of the mass M c . The determination of explicit dependence 
of M s on M c requires the introduction of a specific microscopic model, which is beyond the 
scope of this paper. We point out that, in the static limit u — > 0, and for momenta q <C 1/a, 
the susceptibilities (0) have the form of the Ornstein-Zernicke critical correlation function. 

Close to criticality (T, Mj « t) a perturbative analysis is in principle not justified. 
Nonetheless, to capture the essential features of the single-particle spectra when electrons 
are coupled through a singular effective interaction, we perform our calculations within 
perturbation theory. The agreement with the full calculations |25j in the static limit, as 



well as the claimed smallness of vertex corrections |J, provide a support to the substantial 
validity of our approach. Since a full microscopic derivation of the parameters appearing in 
our model would not be appropriate in this context, we rather follow a phenomenological 
analysis, based on the interpretation of experimental results. 

The lowest order in perturbation theory gives an electron self-energy 

E(k,e) = E c {k,e) + 3Z s (k,e) -S/i, (3) 
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where a factor of 3 in front of S s comes out from the sum over the three spin components 
j = 1,2,3, which are equivalent in the paramagnetic phase. We discuss the single-particle 
spectral properties by means of the spectral density 

A(k e) = 1 |ImE(MI m 

K,) 7r[£-£ fc -Re£(fc,£)] 2 + [Im£(fc,£)] 2 ' K) 

The correction 8[i to the chemical potential is fixed, for a given set of coupling constants 
{Xj}, by the condition 

, / , TTTTv / A(k,e)f(e)de = l-x, (5) 

-w/aJ-TT/a \LTXjay J-oo 

where x is the hole doping per unit cell with respect to half filling, and a factor of 2 in the 
l.h.s. accounts for spin degeneracy. The self-energies appearing in (|3|) are calculated at a 
fixed value of the bare chemical potential fi, which corresponds to the same hole doping x. 
The explicit expression for the imaginary part of the self-energies appearing in ([3]) is 

T v n. n i , ,\ r'° r'° *K -&■][/(&-) 

s) = X^tsgnie) J ^ (2T/a)2 [e _ h , p + sl 2 (k _ kJ) ■ (6) 

where f(e) = [e £ ^ T + is the Fermi function, and b(e) = [e £ ^ T — is the Bose function. 
Within our approach the momentum cut-off q ma x = it/a is provided by the underlying 
lattice. Since, however, we assumed the specific form (Q) for the susceptibilities, which 
is expected to hold for momenta q ~ Qj, we introduce in @ a smooth cut-off function 
C^{k — k! — Qj) = exp(— 7fc_fc'-Q J /A 2 ), whenever a quantitative improvement is needed to 
compare with experimental results. We point out that the presence of a finite cut-off A does 
not change the qualitative behavior of the quantities under discussion. 

At T = the integral over k' in (|6]) is restricted to the region < < min(e, ^ max ) for 
e > and max(e, £ m j„) < < for e < 0, where £ m j n = — 2i — 4i' — /z and £, max = It — At' — fi 
are the minimum and the maximum over k of (in the parameter range — 1 < 2t' /t < 0). 

The real part of each self-energy is obtained as the Kramers- Kronig transformation 
(KKT) of the corresponding imaginary part (|6]), 

p r+oo 

ReS,(fc, e) = — / ImS,(A;, u)sgn(uj) . (7) 

7T J-oo UJ — e 
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Equations @ and (0) are the starting point to obtain single-particle physical properties 
such as the on-shell inverse scattering time, the spectral density ([|) and the characteristic 
features of the FS which will be discussed in the following sections. 



III. ON-SHELL INVERSE SCATTERING TIME. 

The inverse scattering time provides an indication of the violation of the normal FL 
behavior near criticality (i.e. for Mj,T <C t). At the lowest order in perturbation theory 
the on-shell inverse scattering time for a quasi-particle with energy > at T = is 

(1) =-2ImE i (k,£ k ) = 2\ i u i t [ [ , ^f' ? 7 (8) 

\n) 3 A ,qk} 3 3 J J k >: o<e fc ,<^ (27r/a) 2 [& - &f + Oj(k — k') y } 

for each scattering channel j. In this section, as a matter of illustration, we consider the effect 
of a single quasi-critical mode, dropping the index j. We discuss the generic behavior of l/r k 
and individuate the regions where the strongest violations to the FL behavior take place. 
When more than one quasi-critical mode is considered, the resulting behavior is dominated 
by the most singular contribution. For a given k the two terms in the denominator of (|8|) do 
not vanish simultaneously unless = £k-Q- Then the most important contribution to (§) 
comes from a region around the point k' = k — Q, where = and Q(k — k') = M. We 



take n(k — k')~M + va 2 (\k — Q\ 2 (f) 2 + p 2 )/2 where p = \k'\ - \k - Q\ ~ (£ k , - t h - Q )/v k . 



Qi 



\p\ «C \k — Q\, Vk-Q is the velocity of electrons at the point k — Q, and <fi is the (small) 
angle between k' and k — Q. The low-energy behavior of l/r k is most singular when k varies 
along the line = £k-Q, which meets the FS at the "hot-spot" where kp = kus; such that 
£k HS = £,k HS -Q = 0- Then we take Q(k — k') ~ M + z/a 2 |A;^5 — Q\ 2 <p 2 /2 and we perform the 
integral over p in (|5D, finding 

1 2\tis\k HS -Q\a 2 ^/f k 



Tk V k HS -Q J0 



d^log 



x + /M + va 2 \k H s - Q\ 2 G 2 



(9) 



where = (/)/\^ k and the upper limit is extended to oo, to extract the leading behavior. At 
the QCP M = and the integration over 9 yields the non-FL behavior [EH] 
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1 ^ (io) 

When M > two different regimes exist. For ^> M the behavior ((L0|) is again found, 
whereas for ((.<M the mass term prevents the denominator in from vanishing and the 
FL behavior [^6 



1 XtaVuM ( i 



2 



(11) 



T v k HS -Q \ M J 

is recovered. When k approaches kp along a line on which ^ ^ £k-Q, we take — k') ~ 
M* + z/a 2 |fcp — Q\ 2 <f) 2 /2, where M* = M + va 2 £,l F _ Q /vl F _ Q is, in this case, the energy scale 
which separates the FL regime from the anomalous regime and is in general finite at the 
QCP, leading to a FL behavior for < £ K < M*. In this case Eq. ([□]) holds with M -> M*, 
kns ~^ kp- The maximum violation of the FL behavior is found for M* = 0, i.e. M = 
and ^a: f -q = (i.e. k F = k HS ), in which case Eq. (|10|) holds down to = 0. 

We study the dimensionless quantity l/Az/r fc as a function of ^k/imaxi according to 
in the case when k approaches kjjs as ^ — > 0. In this case the asymptotic behavior ([T0|) 
and (|Tl|) are recovered according to the value of the energy scale M (i.e. the distance from 
the QCP). For the sake of definiteness we fix the parameters of the electronic spectrum as 
t'/t = —0.25, and the parameters of the fluctuation spectrum as vjt = 5.0, Q = (ir / a, n / a) . 
For this value of the characteristic wave-vector the hot spots within the first Brillouin zone 
can be determined analytically as the intersections between the FS = 0) and the four 



lines k y = ±n/a±k x . One hot spot is kjis — (a _1 arccosy /i/4t', n/a — a _1 arccosy/i/4t'), and 
the others are found by applying all the symmetry transformations 1Z of the point group. 
We leave the ratio [i/t (or equivalently the hole doping parameter to change the 

position of the Fermi level and the location of the hot spots, and to discuss the resulting 
different regimes. We also let M/t vary in the range [10~ 6 ; 10" 1 ], to explore the cross-over 
from the non-FL regime to the FL regime for k approaching kns- 

In Fig. |](a) we take \ijt = —0.9 (corresponding to a hole doping x = 0.17), i.e. kns — 
(0.32/ a, 2.82/ a). As M/t is reduced (by a factor of ten for each curve from bottom to top), 
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the square-root behavior ([T0|) extends over a wider range of energies. In particular, within 
the energy range considered in Fig. |l|(a), the bottom curve (M/t = 10 _1 ) displays a FL 
behavior (|TTD, whereas the top curve (M/t = 10 -6 ) displays a non-FL behavior (|10D. The 
intermediate curves show the crossover from one regime to the other. 

The above analysis has to be refined if the point kns~Q on the FS is close to a singular 
point of the electronic spectrum so that Vk HS -Q is small. For instance, the characteristic 
wave- vector Q s = (ix/a,it/a) of the AFM fluctuations connects the M points k = (ir/a,0) 
and k' = (0,—ir/a), or equivalently (0,7r/a), of the Brillouin zone, where saddle-point van 
Hove singularities exist at an energy ^vh = At' — fi. A more accurate calculation is needed, 



as the chemical potential approaches /i = At' and Vk F -Q ~ y\fM — fi\, making the estimates 
( [10| , [□]) meaningless. In the limiting case /x = p,, v kp _Q = and we find 



8\vta 2 I dy I dx- — , (12) 



V 2m x ^ 2m y ) ^ 

where l/m x = a 2 (2t + it'), l/m y = a 2 (2t - At'), x = hfj\f^ and y = (k' y - 7r/a)/Vffc. At 
the QCP M = and fll2| ) is independent of indicating an even stronger violation of the 
FL behavior when the hot spot coincides with a saddle point of the electronic spectrum. 
For M > and « M we take x = r cos 9, y = r sin 9 and z = — ^f-) r2 = A e r 2 , 



where Ag ~ 77(6*0 — 9), 9 = ar ct an \Jm y /m x and rj is a suitable constant, so that (12) 
becomes 

— ~ 4Ai/ta 2 / — / dz * — log — . 13 

Tfc io A e i (l_ z )2 + g) 2r / V M / V&/ 

and the van Hove singularity introduces a logarithmic correction to the FL behavior. In Fig. 
|l](b) we take /i/t = p,/t = —1.0 and plot 1/XuTk as a function of £fc/£max, according to 
for approaching fcijs = (vr/a, 0). 

As M/t is reduced (by a factor of ten for each curve from bottom to top), the anomalous 
^-independent behavior extends over a wider range of energies. In particular, within the 
energy range considered in Fig. |l|(b), the bottom curve (M/t = 10 _1 ) displays a FL behavior 
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(fyf), whereas the top curve (M/t = 10~ 6 ) displays a non-FL ^-independent behavior. The 
intermediate curves show the crossover from one regime to the other. 

The enhancement of the phase-space available for scattering makes quasi-particles ill- 
defined in the vicinity of the hot spot. Moreover, as we shall show in the following, quasi- 
critical scattering near the hot spots does not only broaden the quasi-particle peak, but also 
suppresses of the polar behavior of the electron Green function near the Fermi energy. 

In the model considered in this paper, where quasi-critical charge and spin fluctuations 
are simultaneously present, the hot spots corresponding to the different modes co-exist, 
leading to an extension of the anomalous behavior observed along the FS. This fact makes 
less stringent the objection in Ref. |26|] about the limited effect of isolated hot spots on the 
behavior of the normal phase. 



IV. QUASI-PARTICLE SPECTRA 

In this section we discuss the properties of the self-energy (0,0) at T = 0. For the sake 
of clarity, we focus again on a single quasi-critical mode, dropping the index j. It must, 
however, be borne in mind that, once the relevant features are individuated, the physics 
of the system with coexisting charge and spin fluctuations results essentially (though not 
exactly) from a superposition of the separate effects. 

The imaginary part (D), at a fixed k, is characterized by different features depending 
on the value of the energy e compared to characteristic energy scales which we discuss in 
the following. For the sake of the discussion, in Fig. |2|, we fix the parameters as t = 200 
meV, t' = —50 meV and fi = —180 meV (corresponding to a hole doping x = 0.17) for 
the electronic spectrum, M = 2 x 10~ 4 meV, v = 1000 meV and Q = (ir/a,7r/a) for the 
fluctuation spectrum and we take a dimensionless coupling constant A = 0.2. We also fix 
the external momentum k = (1.32/a, 1.32/a), so that ~ —10 meV and £k-Q — 390 meV. 

At energies e > ^ max or e < £ m i„ the integral in (|6]) extended over a domain which 
is independent of e, so that ImE ~ —C/e as soon as |e| ^> max(£ ma2 ., — £ m j n , 2u), where 
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C = ~Xtv[l +x sgn(e)], x being the hole doping per unit cell. Broad dispersionless maxima 
for |ImE| are found at energies e ~ min(£ m j n , —2v) and e ~ max(£ ma;r , 2v). Since in the case 
of Fig. ^(a) 2v > £ max , —£min, the maxima are found at energies e ~ ±2u = ±2000 meV. 

Within the relevant range £, min < e < t; max , two characteristic energy scales exist, M and 
£/c_q. The first controls the low-energy behavior of ImE(&, e —>■ 0), which may be obtained 
by generalizing the methods discussed in the previous section, to the case when e varies 
independently of k. A FL behavior (~ e 2 ) is found for \e\ <C M* = M + va 2 {^k~Q/vk-o) 2 /2 
and a non-FL behavior (~ \J~W\) is found for M* = (i.e. M = and ^-q = 0). The case 
shown in Fig. a corresponds to a low-energy FL behavior . 

A particular role is played by the energy e = £k-Q, since both terms in the denominator 
of (|6|) vanish at k' = k — Q (when M = 0). A peak, or a shoulder partially merged in the 
broad background structure discussed above is present in ImE, depending on the parameters 
of the model. In Fig. 0(a) this point is marked by a diamond at an energy e ~ 400 meV. 
This dispersing structure follows the shadow band ^-q with varying k. As £&_q — > the 
peak at e — £k-Q is suppressed, due to the fact that ImE must vanish at the Fermi energy, 
but when M = the usual FL behavior is turned into an anomalous square-root behavior 
which is present all along the curve £k-Q = 0, usually called shadow FS. However it must 
be pointed out that such an anomalous behavior shows up in the quasi-particle properties 
only when the quasi-particle peak crosses the Fermi energy at the shadow FS, i.e. in the 
vicinity of the hot spots, where ~ £,k-Q — 0. Away from the hot spot, spectral weight is 
transferred from the quasi-particle peak to the shadow peak, leading to an enhancement of 
low-laying spectral weight along the shadow FS. 

The real part of the self-energy, including a correction 5fi ~ —100 meV to the chemical 
potential, found by solving Eq. ([5]) for the present set of parameters, is shown in Fig. ^ b. 
The presence of a peak in ImE at an energy e = £&_q produces a corresponding feature in 
ReE [marked by a diamond at e ~ 400 meV in Fig. 0(b)] which, unlike the case of Ref. JTS 



is not symmetric around e = £k-Q, due to the vanishing of (|5|) at e = 0. As it will become 
clearer in the discussion of the quasi-particle spectra at the end of this section, the resulting 
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suppression of the quasi-polar structure of ReE turns the shadow-bands found in Ref. fI5 



into incoherent resonances at energies e ~ £k-Q- We also point out that, in correspondence 
of the low-energy anomalous behavior ~ v|e|, the suppression of the shadow feature in ReE 



is weaker than for the corresponding FL behavior (~ e 2 ). 

In the regime of FL violation, at low energies ImE = — v4sgn(£:) -y/jej and, by KKT, 
ReE = — Asga(e) \J\e\ + B , where A, B are constants, which depend on the parameters of the 
model. The electron Green function looses then its polar structure. If one however insists in 
assigning a polar structure to the Green function, the wave-function renormalization factor 
Z = {1 — [ReE(e) — ReE(0)]/£:} _1 ~ vanishes at the Fermi energy. The quasi- 



particle lifetime at the hot spots, 1/T qp = 2Z|Im£| = 2|£*.|, turns out to be independent 
of the constants A, B. This behavior is a signature of the importance of considering the 
appropriate renormalizations when studying the transport properties in the presence of 
singular scattering. Thus this problem turns out to be more involved than usually considered 



The single-particle spectral properties are analyzed by means of the spectral density (|]). 
The quasi-particle peak is located, for a given k, at the intercept between the straight line 
y = s — ^k and the curve y = Re£(fc, e). This intersection, in the weak-coupling limit A 1, 
is located at an energy close to [e ~ —40 meV in Fig. 0(b)]. There is, however, a transfer 
of spectral weight to an energy e ~ ^k-Q-, where the denominator in (|J) has a local minimum 
[this point is marked by a diamond in Fig. @(b)]. This incoherent spectral weight follows 
the dispersion of £,k-Q (or of all the ^-q 3 in the general case when more fluctuating modes 
are present) giving rise to dispersing shadow resonances, which substitute in our approach 
the shadow bands found in Refs. fl6|jr5|. As an example, in Fig. |3| spectra are shown along 
the TX direction in k space, away from the hot spots, for the same set of parameters as Fig. 
and momenta kp = (£ir/50a, £tt/50ci), with I = 15, 17, 19, 21. In this region of the Brillouin 
zone both the quasi-particle peak and the shadow peak have comparable intensities. As 
£,k-Q and £fc get closer to each other, an increasing amount of spectral weight is transferred 
from the quasi-particle peak to the shadow peak. 
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V. SPECTRAL DENSITY AND FERMI SURFACE IN ARPES EXPERIMENTS. 



In this section we use the general results obtained above to discuss recent ARPES exper- 
iments Bi2212 near optimal doping p0tpT| , p4tp7| , p8[ . The increasing experimental evidence 



for a stripe phase [|rj|,[nj, reinforces our idea of having a charge modulation which enslaves 
a spin modulation, leading to the model ([I]) with coexisting charge and spin fluctuations. 
Thus, we deal here with the full self-energy (|3]) associated with critical modes of both types. 
We point out again that, when more than one mode is considered, the resulting effects are es- 
sentially a superposition of those discussed in the previous sections for a single generic mode. 
The presence of AFM fluctuations triggered by the incipient stripe instability is described, 
within our model, by the condition that both modes are quasi-critical, i.e. M c < M s <C t. 
We neglect in the following the possibility for a modulation of the AFM characteristic wave- 



vector induced by stripe formation and the related splitting |10| , which adds minor changes 
to the quantities under discussion, and does not appear to be resolved by ARPES data. 
We make a comparison between the spectral density (HI) and the experimental energy 



distribution curves (EDCs) |29| and we analyze the shape and properties of the FS. To 
reproduce typical conditions in ARPES experiments we introduce a convoluted spectral 
density 

A R (k, e) = / de'A(k, e')f(e')S R (e' - e) (14) 

J — oo 

which takes care of the absence of occupied states above the Fermi energy, through the 
Fermi function, and of the experimental energy resolution R, through the resolution function 



£r(£) — exp(— e 2 /2R 2 ) / \J2ttB? . As it is customarily in ARPES experiments, the resolution 
in momentum space is assumed to be finer than the energy resolution, so that no convolution 
over k is performed in (|14]). A comparison is then possible between the experimental EDCs 
and the convoluted density ( |T4T) as a function of e at given k. The self-energy which enters 
the expression (^) is evaluated here with a smooth cut-off A = 0.4, which selects transferred 
momenta close to the characteristic wave-vectors Q c , Q s , to improve the agreement with 
experimental results [see Eq. (H) and the subsequent discussion]. 
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A particular care is required to discuss the experimentally determined FS. Indeed, the 
presence of the shadow peaks and of the associated incoherent spectral weight at low 
energies along the shadow FS, makes the quasi-particle FS [determined by the equation 
£fc + ReS(A;,£: = 0) = 0] not appropriate. The ARPES EDCs reveal the FS via the zero- 
energy crossing of the dispersing intensity peaks, regardless of their coherence. More recently 
the FS was accurately determined as the distribution of spectral weight within an energy 
window around the Fermi energy, by means of the angle scanning photoemission (ASP) ||20|| . 
Once again the experimental distinction between the main and the shadow FS is lost. Thus, 
we associate with each k point the integrated spectral weight of low-laying occupied states 

r w 

p?= A(k,e)f(s)ds, (15) 
J-w 

where W <C t and we compare the distribution of pY with the FS observed in ASP 



experiments ||31 ] . We point out that, in the absence of an anomalous transfer of spectral 



weight to low energies, such a determination of the FS coincides essentially with the theo- 
retical definition in terms of quasi-particle, at least for reasonably small W. 

We want to discuss the evolution of the band structure as determined by following the 
intensity peaks in ARPES experiments and the properties of the FS, which is determined 
with more details in ASP experiments. To reproduce the experimental conditions in ARPES 
experiments |2l| we take T = 10 meV and R = 15 meV in ([14]) . To reproduce the conditions 



in ASP experiments |20] we take T = 25 meV and W = 25 meV in (|i~5|). 



The parameters of the free-electron band will be fixed to reproduce the band structure, 
the FS and the electron filling of Bi2212 at optimal doping, i.e. t = 200 , t' = —50 meV, 
H = —180 meV, corresponding to a hole doping x = 0.17. 

The experimentally observed FS in Bi2212 is approximately mirror-symmetric with re- 
spect to the TX(Y) axes, but not with respect to the TM(Mi) axes ||20|| . An additional 
feature of the experiments is the presence of a dispersing band along the TM direction, 
crossing the Fermi energy at hp = (0, 0.27r/a), which is absent along the TMi direction PS 



Within the stripe-QCP scenario, the symmetry of the FS is expected to be lowered by the 
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interaction between electrons and charge fluctuations, characterized by an incommensurate 
wave-vector Q c , whereas in the presence of spin fluctuations only, the excitation spectrum 
and the FS are expected to have the full symmetry of the lattice, due to the peculiar com- 
mensurability of the characteristic wave- vector of the Neel AFM structure, Q s = (it /a, ir/a). 
Within our model, in order to preserve the symmetry k — > —k, we always perform an average 
over ±Q C , i.e. we take a self-energy E c = |[S C (Q C ) + E c (— Q c )]. Other point symmetries of 



the original system may be broken according to the direction of Q c , which is system [JTO 
and model dependent |33| . 



Based on the symmetry of the FS observed in Bi2212 f20| , for this specific system we 
take Q c ~ (OAir/a, —OAir/a), directed along the ry direction. 

We discuss in some detail the case v c = v s = 200 meV, M c = M s = 10 meV and A c = 1/2, 
X s = 1/6 (to which Fig. prefers). The correction 5fi to the chemical potential, obtained by 
solving Eq. (||) for the present set of parameters, is 5fi = —40 meV. 

The single-particle spectra are characterized by the presence of shadow peaks associated 
with ±Q C and Q s . The two shadow peaks, associated with £,k±Q c , and the corresponding 
branches of the shadow FS (£a,±q c — 0) do not coincide, as it is instead the case for the 
highly commensurate vector Q s (which is equivalent to — Q s ). 

In Fig. |i| the EDCs along the TM direction are reported in the panel on the right. The 
momentum k is increased from top to bottom by a uniform step Ak = 7r/8a. The quasi- 
particle peak moves from the left to the right towards the Fermi energy and interferes with 
the shadow peak associated with +Q C (see the third curve from the top) which moves initially 
to lower energies (fourth curve) and then again towards the Fermi energy. This shadow peak 
appears in Fig. f| as a shoulder merged in the quasi-particle peak in the seventh curve from 
the top. 

The shadow peak associated with Q s is located at high energies above the Fermi level, 
near the T point. However, as k is increased, this peak moves towards the Fermi energy, while 
the quasi-particle looses spectral weight, which is transferred to the shadow peaks. In the 
last three curves in Fig. f| the shadow peak associated with Q s appears as a broad incoherent 
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peak moving towards an energy e ~ —150 meV while approaching the M point. We associate 



this peak with a corresponding feature observed in Bi2212 at an energy e ~ —200 meV jzT 

For a comparison we report, in Figs. [5] and || respectively, the effects due to spin and 
charge modes when considered separately. In the right panel of Fig. [5], from top to bottom, 
the quasi-particle peak looses its weight as it moves towards the Fermi energy, while near 
the M point a broader shadow structure still appears below the Fermi energy. Recently this 
feature was discussed within the NAFL scenario in Ref . f25| , where a full re-summation of the 



diagrammatic expansion in the static limit was performed. We point out, however, that these 
results match those obtained within our simpler perturbative approach. In the right panel 
of Fig. |6| we report the single-particle spectra in the TM direction in the presence of charge 
fluctuations only. The first six curves from the top are very similar to the corresponding 
curves in Fig. U since there the spin-related shadow peak is located at much higher energy 
above the Fermi level. The shadow peak related to +Q C appears to cross the Fermi level at 
kp — (2.4/a, 0) in Fig. |6|. In the last three curves, due to the absence of the effects related 
to spin fluctuations, the quasi-particle peak never reaches the Fermi energy along the 7M 
direction. The van Hove singularity is located at ^vh — ~ 50 meV. 

The description of the additional feature observed only in the TM direction [^8|] would re- 
quire the introduction of an additional charge modulation wit a wave- vector Q' c ~ (0, 27r/3a). 
We do not discuss these finer details here. 

We now discuss the properties of the FS according to the distribution (|15|) of low-laying 
spectral weigh. The main features of the FS, reported in the left panel of Fig. [|, are an 
inhomogeneous distribution of spectral weight, which is maximum near the diagonals of the 
Brillouin zone, and an asymmetric suppression of spectral weight near the M(Mi) points, 
as experimentally observed pD| . The first aspect is related to the evolution of the FS in 
the underdoped regime, as the temperature is lowered down to the superconducting critical 
temperature [0. We point out that, since we are at a fixed doping x = 0.17, and since our 
model (|TJ) is inadequate to discuss the underdoped regime, where the spectrum of charge 
fluctuations is presumably modified by the presence of preformed Cooper pairs, we can only 



describe the onset of this evolution. The second aspect is understood as a co-operative effect 
of spin fluctuations, which lead to a symmetric suppression (see the left panel in Fig. |5|, 
where spin fluctuations only are considered), and charge fluctuations, which introduce an 
asymmetric modulation related to the direction and magnitude of Q c (see the left panel in 
Fig. where charge fluctuations only are considered). 

The branches of the shadow FS due to charge fluctuations, which are clearly visible in 
Fig. H, where they cross transversally the YM{M\) directions in the vicinity of the M(M\) 
points, interfere with the branches due to spin fluctuations, when both mode are considered 
as in Fig. |j. 

Another feature of the FS, related to spin fluctuations, is the possibility for the presence 
of hole pockets around the points (±7r/2a, ±7r/2a) of the Brillouin zone. According to our 
interpretation of the experimentally observed FS we choose a set of parameters which leads to 
very weak hole pockets, which are well below the threshold fixed for the representative points 
in Figs. f| and |5|, and are not reported. Larger \t'/t\, and/or larger u s /t, and/or a sharper 
cut-off A, and/or a smaller energy window W, make this phenomenon more pronounced/, 
Whenever present, the hole pockets are not associated with a topological change of the quasi- 
particle FS, as in Ref. JJ, but rather correspond to a distribution pY which is characterized 
by local maxima at the FS (£& = 0) and at the branch of the shadow FS (£,k~Q s = 0), 
and minima at (±7r/2a, ±7r/2a). The experimental situation is still controversial within this 
respect. Shadow bands shifted by a wave- vector Q s and/or the related hole pockets were 
reported in Refs. [20,|24J35| but were not detected in Refs. pl|,|34| . 



VI. CONCLUSIONS 

We studied the changes in the single-particle properties due to the coupling of electrons 
to quasi-critical charge and spin fluctuations in the quantum critical region around a stripe 
QCP. Violations to the normal FL behavior appear in the on-shell inverse scattering time, at 
those points of the FS connected by a characteristic wave-vector of the critical fluctuations 
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(hot spots). The violations are stronger when a hot spot is located near a singular point 
of the electronic spectrum, leading to a finite inverse scattering time. This result may be 
relevant for the physics of high-T c cuprates, where extended van Hove singularities have 



been repeatedly observed close to the Fermi energy at the M points, in the presence 
of critical fluctuations characterized by either a large (~ Q s ) or a small Q c wave-vector, 
leading to hot spots in the vicinity of the M points. 

The study of the quasi-particle spectra showed that the violation to the FL behavior are 
associated with the transfer of spectral weight from the quasi-particle band to incoherent 
shadow peaks. These feature are dispersing and show up in the ARPES EDCs, though they 
do not correspond to poles of the electron Green function. As the shadow peaks approach 
the Fermi level, they produce an enhancement of low-laying spectral weight which is seen 
shadow FS within a finite energy resolution. 

The evolution of the FS is thus associated with the change in the distribution of the 
low-laying spectral weight and not with the topological modification of the quasi-particle 
FS proposed in Ref. |J. 

The quasi-particle spectra and the FS obtained within the model with coexisting charge 
and spin fluctuations capture the relevant features observed in ARPES experiments on 
optimally doped Bi2212, and namely the presence of broad features in the EDCs (besides 
the quasi-particle peak) the inhomogeneous distribution of spectral weigh along the FS and 
the asymmetric suppression of spectral weight around the M(M X ) points of the Brillouin 

0- 



zone 



On the other hand we suggest that the ARPES experiments may be used to extract 
informations about the parameters entering in the fluctuation spectra (Q), which can then 
be tested consistently within the theory. In particular, the direction and magnitude of Q c 
might be deduced from the modulation in the suppression of spectral weight along branches 



of the FS, as suggested in Ref. po| , or from a systematic, and still lacking, investigation of 
the related shadow peaks. 
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FIGURE CAPTIONS 

Fig. |1|: Dimensionless on-shell inverse scattering time at the hot spot, as a function 
of ik/imax- (a) The parameters of the electronic spectrum are t'/t = —0.25, fi/t = —0.9 
(corresponding to a hole doping x = 0.17); the parameters of the fluctuation spectrum are 
v/t = 5.0, Q = (vr/a, vr/a), while M/t varies in the range [10~ 6 ; 10 _1 ] increasing by a factor 
of 10 for each curve from top to bottom. The hot spot is located at kns — (0.32/a, 2.82/a). 
(b) Same parameters as in (a), except for the chemical potential n/t= —1.0 (corresponding 
to a hole doping x = 0.22). The van Hove singularity in the electronic spectrum is located 
at the Fermi level. The hot spot is located at kns — (7r/o> 0). 

Fig. |2|: Self-energy as a function of the external energy e. The parameters are the same as 
in Fig. [l|(a), with t = 200 meV, to fix the energy scale, and M/t = 10~ 6 . The dimensionless 
coupling is A = 0.2 and the momentum is k = (1.32/a, 1.32/a). For these values of the 
parameters ~ —10 meV and £k-Q — 390 meV (marked by a diamond on the energy axis), 
(a) Imaginary part. The broad dispersionless maxima are located at e ~ ±2000 meV. The 
shadow feature is located at e = £fc-Q- (b) Real part, including a correction to the chemical 
potential S/j, ~ —100 meV. The shadow feature is located at e = £,k-Q- The quasi-particle 
peak corresponds to the intercept between the solid line e — and the curve ReS(fc, e), and 
is located at an energy e ~ —50 meV. The dashed lines represent the axes. 

Fig. 0: Spectral density for kt = (in/ 50a, in /50a), with £ = 15,17,19,21 (dot-dashed, 
dotted, dashed and solid lines respectively). The other parameters are the same as in Fig. 
0. The quasi-particle peak is moving from the left to the right (towards the Fermi energy) 
as I is increased. The shadow peak is moving from the right to the left as I is increased. The 
position of the shadow band, e = Cfc-Q, is marked by diamonds on the energy axis. Spectral 
weight is transferred from the quasi-particle peak (on the left) to the shadow peak (on the 
right) as the two peak get closer to each other. 

Fig. |j: LEFT: Distribution of low-laying spectral weight pY within the Brillouin zone 
in the case of electrons coupled to both charge and spin fluctuations. The weight is reduced 
by a factor of 2 as the size of the black squares is reduced. Four classes are shown. The 
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weight less than 1/16 of the maximum is not reported. The parameters of the electronic 
spectrum are t = 200 meV, t' = 50 meV, fx = —180 meV (corresponding to a hole doping 
x = 0.17). The parameters of the fluctuation spectra are u S}C = 200 meV, M S}C = 10 meV, 
Q s = (vr/a, 7r/o) and Q c = (OArc/a, — 0.47r/a). The dimensionless coupling constants are 
A s = 1/6, A c = 1/2 and the dimensionless cut-off is A = 0.4. The correction to the chemical 
potential is 5 fx ~ —40 meV. The integrated spectral weight pjf is calculated according to 
(PI) with W = 25 meV and T = 25 meV. RIGHT: Convoluted spectral density A R (k,e), 
calculated according to ([14]) with R = 10 meV and T = 10 meV, along the TM direction 
for the same parameters as in the left panel. The momentum is uniformly increased from 
top to bottom. 

Fig. |5|: LEFT: Distribution of low-laying spectral weight pjf within the Brillouin zone 
in the case of electrons coupled to spin fluctuations only. The weight is reduced by a factor 
of 2 as the size of the black squares is reduced. Four classes are shown. The weight less 
than 1/16 of the maximum is not reported. The parameters of the electronic spectrum are 
the same as in Fig. The parameters of the spin-fluctuation spectrum are v s = 200 meV, 
M s = 10 meV, Q s = (7r/a, 7r/a). The dimensionless coupling constant is A s = 1/3 and the 
dimensionless cut-off is A = 0.4. The correction to the chemical potential is Sfi = —70 meV. 
The integrated spectral weight is calculated according to dig) with W = 25 meV and 
T = 25 meV. RIGHT: Convoluted spectral density A^(k,e), calculated according to ([T4|) 
with R = 10 meV and T = 10 meV, along the TM direction for the same parameters as in 
the left panel. The momentum is uniformly increased from top to bottom. 

Fig. H LEFT: Distribution of low-laying spectral weight pjf within the Brillouin zone 
in the case of electrons coupled to charge fluctuations only. The weight is reduced by a 
factor of 2 as the size of the black squares is reduced. Four classes are shown. The weight 
less than 1/16 of the maximum is not reported. The parameters of the electronic spectrum 
are the same as in Fig. [|. The parameters of the charge-fluctuation spectrum are v c = 200 
meV, M c = 10 meV, Q c = (0.47r/a, — 0.47r/a). The dimensionless coupling constant is 
A c = 1/2 and the dimensionless cut-off is A = 0.4. The correction to the chemical potential 
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is negligible. The integrated spectral weight pY is calculated according to ( [L5|) with W = 25 
meV and T = 25 meV. RIGHT: Convoluted spectral density A R (k,e), calculated according 
to (|T4| ) with R = 10 meV and T = 10 meV, along the TM direction for the same parameters 
as in the left panel. The momentum is uniformly increased from top to bottom. 
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FIGURES 




FIG. 1. S. Caprara et al. - Single-particle properties of a model ... 



29 




FIG. 2. S. Caprara et al. - Single-particle properties of a model ... 
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FIG. 3. S. Caprara et al. - Single-particle properties of a model 
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FIG. 4. S. Caprara et al. - Single-particle properties of a model ... 
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FIG. 6. S. Caprara et al. - Single-particle properties of a model ... 
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